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1.0  INTRODUCTION 

This  report  documents  the  findings  of  an  experiment  which  placed  a  wide¬ 
band  (40  to  4000  Hz)  set  of  vertical  acoustic  antennas  in  The  Tongue  of  the 
Ocean  (TOTO)  at  a  site  adjacent  to  the  Andros  Island  eastern  escarpment.  This 
report  details  the  instrumentation  and  data  processing  procedures  and  docu¬ 
ments  the  findings  of  a  general  analysis.  The  purpose  of  the  experiment  was 
the  initial  field  test  of  a  deployable  acoustic  monitoring  system  (DAMS)  capa¬ 
ble  of  autonomous  remote  operation.  The  success  of  the  experiment  made  avail¬ 
able  a  data  set  that  allowed  the  examination  of  the  vertical  space -time  sec¬ 
ond-order  statistics  of  the  underwater  acoustic  ambient  in  an  area  of  interest 
to  the  Navy.  The  range  of  wind  speeds  (2  to  16  m/s)  was  large  enough  during 
the  experiment  period  (September-October  1988)  to  draw  statistical  conclusions 
relative  to  the  acoustic  ambient  of  the  area  with  particular  emphasis  on  the 
"sea  surface  sound"  that  dominates  the  data.  The  25  data  sets  examined  repre¬ 
sent  approximately  semidiurnal  samples  of  the  environment.  The  primary  objec¬ 
tives  of  this  report  are  (1)  to  document  the  wind  speed  dependence  of  the 
sound  pressure  level  autospectral  density  and  the  acoustic  vertical  direc¬ 
tional  spectra,  (2)  to  identify  the  lower  bounds  for  the  anisotropic  gains 
that  would  be  realized  by  vertical  acoustic  antennas  operating  in  the  area, 
and  (3)  to  report  on  the  industry- caused  acoustic  radiation  monitored  during 
the  experiment.  A  secondary  objective  is  to  document  some  system  aspects  of 
the  initial  deployment  of  the  DAMS. 

The  next  section  of  the  report  provides  all  of  the  fundamental  hardware, 
software,  and  environmental  data  required  to  interpret  the  results.  That  is 
followed  by  a  section  that  notes  some  specific  system  self-noise  limitations 
to  the  end  product  noted  on  the  first  deployment  of  the  system.  The  fourth 
section  covers  the  results  of  the  measurement,  and  the  final  section  summa¬ 
rizes  the  entire  project. 
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2.0  EXPERIMENT  DESCRIPTION 


2.1  IN-WATER  HARDWARE 

A  figure  describing  the  deployable  acoustic  monitoring  system  is  given  in 
figure  1  which  illustrates  the  components  of  the  DAMS.  The  heart  of  the  sys¬ 
tem  is  the  acoustic  module;  it  houses  the  seven  acoustic  antennas  that  cover  a 
frequency  range  of  40  to  4000  Hz.  Directly  above  the  acoustic  module  is  an 
engineering  module  that  houses  signal  processing  units  and  two  array  tilt  sen¬ 
sors.  Above  and  below  these  two  units  are  vibration  isolation  modules  (VIMs) 
that  isolate  the  acoustic  module  from  longitudinal  vibrations  induced  into  the 
system.  A  30-m  double -armored  steel  cable  ; jove  the  upper  VIM  serves  as  an 
electrical  channel  for  the  acoustic  signals,  array  power,  and  array  tilt  sig¬ 
nals  to  the  instrumented  pressure  vessel  (IPV)  which  is  inside  a  syntactic 
foam  collar  that  keeps  the  entire  system  oriented  vertically.  A  depth  sensor 
is  located  in  the  IPV.  The  cable  below  the  system  is  connected  to  a  variable - 
depth,  bottom-moored  system. 

The  acoustic  sensor  system  consists  of  seven  octavely  coalesced  four- 
wavelength  acoustic  line  antennas  covering  the  frequency  range  of  40  to  4000 
Hz.  The  14-bit  hydrophone  signals  from  each  antenna  section  are  transmitted, 
one  antenna  at  a  time,  to  a  subsurface  data  logger  (IPV)  that  synchronously 
records  them  digitally  on  video  cassette  recorders  (VCRs)  along  with  array 
tilt  and  depth  information  as  well  as  the  date  and  time.  A  system  controller 
is  programmed  for  a  wide  selection  of  on/off  times.  There  is  enough  VCR  tape 
and  battery  life  for  a  60-day  deployment  of  alternate -day  samples. 

Acoustic  antenna  details  pertinent  to  this  report  are  given  in  table  1. 

In  all  cases  the  apertures  are  four  wavelengths  wide  at  the  design  frequency. 

The  array  is  bottom  moored  and  maintained  at  nearly  a  constant  depth. 

Two  depths  were  used  in  the  measurement  period:  156  m  and  140  m,  as  measured 
at  the  IPV.  Currents  would  tilt  the  system  and  thus  increase  this  depth.  The 
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INSTRUMENTED  PRESSURE  VESSEL 
(Subsurface  Data  Logger) 

-  CONTROL  ELECTRONICS 

-  DIGITAL  TAPE  RECORDERS 

-  BATTERIES 

-  FLOTATION 

-  DEPTH  SENSOR 


- STEEL  ARMOR  CABLE 


Figure  1. 


Components  of  the  Deployable  Acoustic  Monitor  System 
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Table  1.  Acoustic  Antenna  Details 


Acoustic 

Array 

Section 

Design 

Frequency 

(Hz) 

Line 

Array 

Length 

(m) 

Hydrophone 

Uniform 

Separation 

(m) 

Number 

of 

Hydrophones 

1 

5624* 

0.93 

0.133 

8 

2 

2812 

1.8' 

0.267 

8 

3 

1406 

3.73 

0.533 

8 

4 

703 

7.47 

1.067 

8 

5 

351 

14.94 

2.134 

8 

6 

176 

29.87 

4.267 

8 

7 

88 

59.74 

8.534 

8 

*Note  that 

the  data  recording 

system  limited 

the  upper  frequency 

to  4000  Hz . 

typical  ranges  of  depths  were  156  to  159  m  and  139  to  147  m.  The  center  of 
the  acoustic  aperture  was  60  in  deeper  than  these  depths.  The  array  tilt 
frequency  of  occurrence,  i.e.,  histogram,  is  presented  in  figure  2.  The  fig¬ 
ure  illustrates  that  while  tilt  angles  were  as  high  as  8  degrees  the  much  more 
probable  occurrence  was  for  the  array  to  be  within  3  degrees  of  vertical.  A 
vertically  uniform  current  can  be  inferred  from  these  results  using  the  meth¬ 
ods  developed  for  the  hydrodynamics  of  cables.  For  example,  see  reference  1. 
The  relation  between  current  and  array  tilt  is  also  shown  in  figure  2,  indi¬ 
cating  that  the  current  rarely  exceeded  0.5  knots. 

2.2  DATA  PROCESSING 

The  data  processing  system  is  shown  as  a  data  flow  diagram  in  figure  3. 
The  typical  3-minute,  eight-hydrophone  (32  Mbyte)  space-time  series  coming 
from  the  VCR  tape  playback  in  real  time  is  acquired  by  an  array  processor  in 
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FREQUENCY  OF  OCCURRENCE  OF  ARRAY  TILT 


CURRENT  (knots) 


Figure  2.  Array  Tilt  Measurements 


TR  8579 


Figure  3.  Data  Processing  Flow  Diagram 
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"real"  time.  The  individual  time  series  resident  in  the  array  processor  are 
checked  for  nonstationarities ,  i.e.,  trends  and  "wild"  points,  and  edited  as 
appropriate.  The  cross-spectral  density  function  is  estimated  for  each  pos¬ 
sible  hydrophone  pair.  We  refer  to  this  ensemble  of  cross  -  spectral  density 
functions  as  the  cross-spectral  density  (csd)  matrix.  The  csd  matrix  is  the 
measurement's  fundamental  statistical  description.  The  data  acquisition  sys¬ 
tem  records  all  time  series  using  a  9790.2  Hz  sampling  frequency,  independent 
of  the  acoustic  aperture.  The  spectral  analysis  system  digitally  low-pass 
filters  and  decimates  each  of  the  space -time  series  according  to  its  aperture 
so  as  to  optimize  the  spectral  resolution  for  the  fixed-point  (1024)  fast 
Fourier  transform  (FFT)  size.  Table  2  lists  the  filter  cutoff  frequency  (at  3 
dB)  and  the  frequency  resolution  for  each  acoustic  antenna.  The  spectral  res¬ 
olution  accounts  for  the  1.63  spreading  factor  associated  with  the  "Hann"  tem¬ 
poral  window  used  with  all  the  data.  The  cross-spectral  density  function 
between  the  nth  and  mth  hydrophone  is  estimated  by  using  the  relation 


S 

p 


(*_•")  -  r  j, 


X  (w)X  (w) 

nj  mj 


(1) 


where  Xnj(w)  and  XmJ(w)  are  the  j th  realization  of  the  Fourier  transforms  of 
the  nth  and  mth  hydrophones,  w  is  the  angular  frequency,  M  is  the  number  of 
transforms  averaged,  and  ^  is  the  spatial  separation  between  the  nth  and  mth 
hydrophones.  Only  nonoverlapping  transforms  are  used.  The  actual  number  of 
FFTs  used  in  each  average  varies  with  the  acoustic  aperture,  the  length  of  the 
time  series,  and  the  VCR  tape  operations.  The  typical  number  of  transforms 
for  the  3 -minute  recording  time  for  each  of  the  apertures  is  listed  in  table 
3.  This  table  also  summarizes  the  95%  confidence  interval  for  the  spectral 
estimates  in  each  of  the  acoustic  apertures.  From  a  statistical  sampling  per¬ 
spective  the  estimates  are  quite  good. 
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Table  2.  Spectral  Analysis  Details 


Acoustic 

Aperture 

Number 

Design 

Frequency 

(Hz) 

Cutoff 

Frequency 

(Hz) 

Frequency 

Resolution 

(Hz) 

1 

5624 

4000* 

9.6 

2 

2812 

3944 

9.6 

3 

1406 

1972 

4.8 

4 

703 

1000 

2.4 

5 

352 

500 

1.2 

6 

176 

250 

0.6 

7 

88 

140 

0.3 

*Limitation 

imposed  by 

data  recording 

system. 

Table  3.  Spectral  Analysis  95%  Confidence  Interval  Summary 
(based  on  a  typical  3 -minute  data  acquire) 
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After  the  csd  matrix  is  calculated,  the  spatial  homogeneity  of  the 
measurements  is  verified  by  checking  that  the  csd  values  along  any  diagonal  of 
the  matrix  are  not  different  from  each  other  in  a  statistical  sense.  The 
statistical  uniformity  of  the  matrix  diagonals  is  a  result  of  the  acoustic 
antennas  having  uniformly  spaced  hydrophones  in  a  homogeneous  pressure  field. 

A  matrix  with  this  property  is  termed  a  Toeplitz  matrix.  All  csd  matrices 
calculated  from  a  real  field  variable  have  complex  symmetry  about  the  major 
diagonal.  A  matrix  having  this  symmetry  is  called  a  Hermetian  matrix.  Occa¬ 
sionally  in  the  experiment  the  Toeplitz  property  was  not  found  to  be  true  in 
aperture  7 ,  the  low  frequency  aperture ,  and  in  those  cases  the  offending 
hydrophones  were  omitted  from  the  analysis.  This  point  will  be  discussed  in 
more  detail  in  the  next  section.  Missing  elements  of  the  matrix  caused  by 
omitted  hydrophone  data  are  filled  using  the  Toeplitz  property  when  the 
offending  hydrophone  is  "within"  the  matrix.  If  the  hydrophone  is  an  "end" 
hydrophone,  calculations  are  made  with  the  aperture  shortened  by  a  single 
hydrophone  spacing  rather  than  by  attempting  extrapolation.  Throughout  the 
experiment  the  seventh  hydrophone  (hydrophone  31)  in  the  lowest  frequency 
antenna  (aperture  7)  was  nonfunctioning  and  its  csd  locations  were  filled  as 
discussed  above.  This  process  is  identified  in  figure  3  as  "matrix  editing." 

The  next  step  in  figure  3  is  the  calculation  of  a  wavenumber -frequency 
spectra  from  the  measured  csd  matrix.  The  wavenumber -frequency  spectra  is 
estimated  using  a  modified  form  of  the  Blackman-Tukey  algorithm,2  which 
implements  the  wavenumber  spectra  by  taking  a  Fourier  transform  of  the  spatial 
correlation  function  present  in  the  csd  matrix.  The  modification  accounts  for 
the  complex -number  nature  of  the  csd  matrix  elements  in  place  of  the  real- 
valued  time  series  samples  of  reference  2.  The  algorithm,  which  uses  the 
Hermetian  symmetry  of  the  csd  matrix  and  employs  Hamming  spatial  weighting,  is 


S  (k  ,w)  -  0.23  V(k  ,w)  +  0.54  V(k  ,»)  +  0.23  V(k  ,«) 

p  q  q-X  q  q+1 


(2) 
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where 


V(k  ,w)  -  d 
q 


N-2 


S  (o,w)  +  2  S  Re 

p  r-l 


S  (rd,w)  cos 
p 


rqx 


r 

1 

+  Re 

S2[(N  -  l)d,w 

cos (rx) 

l 

- 

k  -  —  q  (wavenumber) , 
q  Nd 


d  is  the  interhydrophone  spacing,  and 


2 

S  (  )  is  the  diagonally  averaged  csd  element, 
p 

Note  that  the  wavenumber  is  the  line  array  longitudinal  wavenumber. 
Results  of  these  calculations  are  discussed  in  the  next  section. 

The  next  step  in  figure  3  is  the  calculation  of  the  direction  of  arrival 
spectra  (DOAS),  which  is  the  standard  beamformer  output  autospectral  density. 
While  the  csd  matrix  calculation  is  the  fundamental  statistical  quantity  of 
the  experiment,  the  DOAS  is  the  fundamental  observation  of  the  experiment 
because  there  is  no  illuminating  way  to  view  a  csd  matrix.  The  DOAS  on  the 
other  hand  is  a  familiar  and  intuitive  quantity  in  the  study  of  underwater 
acoustics.  It  is  calculated  by  finding  the  weighted  quadratic  form  of  the  csd 
matrix  as  follows 


S 

b 


(3) 
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where 


-M 


w 


wz  exp  j^i^d  cos<f>  j 


w„  exp  ic( 


exp|ic(N  -  l)d  cos^ 

N  [_  S 


S  («)  is  the  NxN  csd  matrix, 
~p 


<t>  is  the  "steered"  elevation  angle  with  zero  degrees  being  straight  up,  and 

w  is  the  real  spatial  weight  applied  to  the  nth  hydrophone. 

N 


A  Taylor  pattern  function  was  designed  for  the  experiment.3,4  The 
spatial  weights  for  the  eight -element  hydrophone  array  are  given  in  table  4. 
The  beam  pattern  functions  are  shown  in  figures  4  and  5  for  antenna  nondimen- 
sional  apertures  of  2  and  4,  with  the  latter  being  the  "design"  frequency. 

The  half-power  beamwidths  for  the  steering  angles  in  the  figures  are  given  in 
table  5.  Note  from  the  figures  that  the  Taylor  pattern  has  nearly  equal  side 
lobes  19  dB  below  the  mainlobe. 
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Table  4.  Taylor  Spatial  Weights 


Hydrophone  Number 

Spatial  Weight 

1 

0.7274 

2 

0.9086 

3 

1.1147 

4 

1.2493 

5 

1.2493 

6 

1.1147 

7 

0.9086 

8 

0.7274 

Table  5.  Taylor  Pattern  Half- Power  Beamwidths 


Steering  Direction 
(Degrees) 


Nondimens ional  Aperture 
(Degrees) 


it  2  •* 


«4" 


0 

83.8 

58.0 

45 

46.0 

20.5 

90 

28.7 

14.2 
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The  next  step  in  the  data  flow  of  figure  3  is  the  formation  of  propor¬ 
tional  bandwidth  spectra.  There  are  two  reasons  for  wanting  to  form  propor¬ 
tional  bandwidths .  First  is  a  reduction  in  data.  Second  is  an  increase  in 
the  statistical  stability  of  the  csd  estimates.  The  first  reason  is  the  more 
compelling  one.  There  are  over  2xl05  complex  words  resulting  from  a  single 
acquire  of  all  seven  antennas.  Much  of  the  data  is  redundant  because  the 
physical  mechanisms  governing  the  process  do  not  change  rapidly  with  fre¬ 
quency.  The  relatively  high-resolution  constant-bandwidth  analysis  over¬ 
samples  the  process  when  this  is  the  case.  Thus  the  csd  matrix  can  be  aver¬ 
aged  over  blocks  of  frequency  with  large  reductions  in  the  amount  of  data 
storage  required  and  with  an  improvement  in  the  quality  of  the  statistical 
estimates.  Because  the  data  need  to  be  archived  and  sorted  with  respect  to 
environmental  parameters,  principally  wind  speed,  there  is  significant  motiva 
tion  to  reduce  the  amount  of  data  that  must  be  sorted.  Thus  we  formed  propor 
tional -bandwidth  csd  matrices.  The  bandwidth  choice  was  arbitrary  because 
these  bandwidths  are  synthesized  from  the  raw  data.  One -third- octave  band- 
widths  were  synthesized  for  both  historical  reasons,  i.e.,  the  prevalence  of 
one- third-octave  analysis  in  the  acoustic  industry,  and  because  the  one-third 
octave  analysis  still  appears  to  over-sample  the  process  we  are  measuring. 

The  next  step  in  figure  3  is  the  eigenstructure  decomposition  of  the  csd 
matrix.5  This  analysis  was  undertaken  with  two  objectives  in  mind.  First, 
this  approach  readily  leads  to  the  identification  of  plane-wave  interference 
in  the  ambient  data.6  The  next  section  of  this  report  discusses  an  extension 
to  this  whereby  the  technique  is  used  to  remove  a  contaminating  interference 
source.  Second,  the  results  of  the  decomposition  are  useful  in  interpreting 
certain  aspects  of  the  data. 

The  "solution"  to  the  eigenstructure  problem  is  well  known7  and  requires 
solving  the  linear  algebraic  equation 

S  (w)H(w)  -  H(u)A(«)  ,  (4) 

-p  -  -  ~ 
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where 


H(w)  is  the  lxN  matrix  of  column  (N- element)  eigenvectors  of  S  (w)  ,  and 

~  ~p 


A(w)  is  the  NxN  diagonal  matrix  of  eigenvalues  [A(w>]  of  S  (w)  . 

“  p 

The  reason  for  finding  the  eigenvectors  and  eigenvalues  is  that  they  represent 
a  decomposition  of  the  csd  matrix  into  a  complete  set  of  linearly  independent 
components.8  The  decomposition  is  explicitly  expressed  as 

N  ^'T 

S  (w)  -  Z  A  (w)H  H  (w)  ,  (5) 

— p  n“l  n  ~-n  — n 


where  An  is  the  nth  scalar  eigenvalue  and  H^w)  is  the  nth  N-element 
eigenvector.  It  will  be  instructive  to  take  the  DOAS  of  both  sides  of  equa¬ 
tion  (5)  by  applying  equation  (3).  This  leads  to  a  decomposition  of  equation 
(3)  into  "eigenbeams" ,  which  are  defined  as 


E  (w) 

n 


A 

n 


-'-Ti  ~n  —  ^  sj 

Q 


(6) 


The  final  step  shown  in  figure  3  is  the  calculation  of  the  vertical 
directional  spectra  from  the  measured  csd  matrix.  Much  has  been  written  about 
estimating  the  directivity  spectra  from  measurements  made  on  linear  hydrophone 
arrays.  The  difficulty  is  that  both  the  beamformer  outputs  and  the  cross- 
spectral  densities  are  related  to  the  directivity  spectra  by  an  integral  equa¬ 
tion  in  which  the  directivity  spectra  is  in  the  integrand.  Inverting  this 
Integral  equation  has  basic  uniqueness  difficulties  when  using  beamformer  out¬ 
put  autospectral  density.*-10  In  the  case  of  making  the  estimate  directly 
from  the  csd  matrix,  an  assumption  that  the  field  is  composed  of  statistically 
independent  plane  waves  is  necessary . u- 12  This  assumption  may  not  be  a  good 
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one  when  significant  energy  is  being  reflected  from  the  ocean  bottom.  Deep- 
ocean  application  of  these  techniques  has  been  made . 11,13,14  A  first-order 
approximation  to  the  directivity  function  can  be  obtained  by  normalizing  each 
beamformer  output  by  its  solid-angle  coverage.6,15,16  More  specifically,  the 
DOAS  may  be  expressed  in  terms  of  the  beamformer  pattern  function  and  the 
desired  directional  spectra  as 


rtf 


2  tf 


d<j>  sin  4>  B' 


:  N' 


(*> 


(7) 


where  B2(^;^s  is  the  beamformer  pattern  function,  assuming  azimuthal  symmetry, 
and  N2(<£)  is  the  vertical  directional  spectra. 

If  the  directional  spectra  varies  slowly  relative  to  the  beam  pattern 
mainlobe  width,  equation  (7)  can  be  approximated  by 

Sb[«>;0s]  -  b2(^s]n2[<*s]  a*(*J  ,  <8> 


where  A V>(<£s)  is  the  solid- angle  coverage  of  the  beamformer  mainlobe,  which  is 
evaluated  using 


AV> 


r  •  + 


2jt 


d<£  sin  <f> 


(9) 


where  and  <f>+  are  the  half-power  angles  associated  with  steering  angle  <j>3. 
The  first-order  approximation  is  realized  by  solving  for  N2($)  in  equation 
(8). 
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This  first-order  approach  yields  reasonably  accurate  results  for  most 
values  of  the  elevation  angle  when  the  acoustic  field  results  from  a  spatially 
diffuse  source  structure.  The  exception  to  this  is  near  the  horizontal  direc¬ 
tion  where  significant  curvature  is  possible. 

2 . 3  EXPERIMENT  ENVIRONMENT 

Figure  6  illustrates  the  measurement  location  relative  to  the  adjacent 
escarpment,  including  the  local  bathymetry.17  The  measurement  site  is  on  the 
eastern  shore  of  Andros  Island  in  The  Bahamas  just  south  of  Fresh  Creek,  as 
shown  in  figure  7.  The  depth  at  the  measurement  site  is  700  meters  (383  fath¬ 
oms),  which  is  on  the  sloping  bottom  shown  in  figure  7. 

The  experiment  took  place  in  September  and  October  1988.  The  near¬ 
surface  sound  velocity -depth  profiles,  typical  of  that  location  and  season, 
are  shown  in  figure  8.  The  two  profiles  were  measured  coincident  with  acous¬ 
tic  measurements.  The  important  acoustic  characteristics  of  the  profiles  are 
(1)  the  sound  velocity  at  the  measurement  depth  is  always  less  than  the  sur¬ 
face  velocity,  and  (2)  the  measurement  depth  is  always  below  the  seasonal 
thermocline  depth,  which  varied  from  30  to  73  meters.  There  are  two 
consequences  to  the  first  characteristic.  First,  near- surface  sound  sources 
contribute  no  energy  to  the  directional  spectra  at  angles  of  less  than  9 
degrees  from  the  horizontal.  Second,  the  near-surface  sound  sources  that  con¬ 
tribute  to  the  measured  sound  field  are  limited  to  a  circular  area  with  an 
approximate  radius  of  1300  meters.  This  limited  area  of  influence  is  shown  in 
figure  6.  The  second  characteristic  indicates  that  the  measurements  never 
include  "surface  duct"  effects. 
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3.0  SYSTEM  SELF-NOISE 


This  report  documents  the  initial  deployment  of  the  DAMS.  NUSC  experi¬ 
ence  with  such  vertical  hydrophone  array  structures  is  limited.  The  purpose 
of  this  section  is  to  document  structural  and  electronic  problems  found  in 
this  deployment.  The  analysis  of  the  data  indicated  three  sources  of  system 
self -noise,  i.e.,  hydrophone  response  due  to  nonacoustic  signals  caused  by  the 
measurement  system.  These  three  sources  of  contamination  are  (1)  common-mode 
electronic  nci.se  generated  by  the  IPV,  (2)  hydroacoustic  noise  caused  by  the 
antenna  structure  interacting  with  ocean  currents,  and  (3)  preamplifier  elec¬ 
tronic  noise. 


The  term  "common-mode  noise"  is  used  to  describe  a  system  electronic 
noise  that  occurs  in  all  hydrophone  channels  with  no  difference  in  phase. 
When  such  a  noise  enters  a  beamforming  algorithm  it  adds  coherently  in  the 


"broadside"  steered  direction,  which  in  this  case  is  an  elevation  angle  of 
90°.  Figure  9,  which  shows  the  direction-of-arrival  spectra  for  the  quietest 
acoustic  ambient  condition  observed,  illustrates  the  common-mode  noise  prob¬ 
lem.  The  2  m/s  wind  speed  condition  shown  in  figure  9  contains  the  lowest 
acoustic  levels  of  the  experiment  and  thus  best  illustrates  the  self-noise 
problem.  At  higher  wind  speeds  the  acoustic  signal  increases  and  partially 
masks  the  self-noise.  Figure  9  depicts,  in  an  elevation  angle  -  frequency 
space,  the  portion  of  the  direction-of-arrival  spectra  contaminated  by  the 
noise.  Figure  10,  which  shows  the  autospectral  density  of  the  beamformer  out¬ 
put  steered  to  an  elevation  angle  of  90®,  gives  another  view  of  the  problem. 

The  common-mode  noise  problem  was  found  to  originate  in  the  instrumented 
pressure  vessel  (IPV)  electronics.  The  battery  supply  voltage  was  used  as  an 
input  to  a  dc-to-dc  converter  and  voltage  regulation  subsystem.  One  of  the 
outputs  of  the  voltage  regulation  system  was  a  12-volt  line  that  drove  the 
VCRs  as  well  as  the  the  hydrophone  array  and  matrix  switch  electronics.  The 
recorders  periodically  "pulled  down"  this  supply  voltage.  This  periodic  sig¬ 
nal,  or  more  exactly  a  harmonic  of  that  signal  which  occurred  in  the  passband 
of  the  system,  appeared  on  the  hydrophone  output  signals.  This  problem  has 
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since  been  corrected  by  including  a  second  dc-to-dc  converter  in  the  IPV  to 
service  just  the  VCRs. 

It  is  clear  from  the  data  that  the  structure  was  undergoing  forced  vibra¬ 
tion  at  certain  times.  Those  vibrations  manifested  themselves  as  very  low 
frequency  signals.  Typically  those  signals  were  limited  to  the  two  end 
hydrophones,  i.e.,  hydrophones  number  1  and  32,  which  are  46  cm  from  the 
acoustic  module  end  terminations.  Those  low  frequency  signals  normally  were 
not  observed  on  hydrophones  2  through  30,  which  are  hydrophones  removed  from 
the  module  end  terminations.  (Hydrophone  31  was  nonfunctional.)  Each  of 
those  hydrophones  is  in  excess  of  9  meters  from  a  module  end  termination.  It 
is  known  that  this  type  of  structure  has  a  resonant  axisymmetric  transverse 
wave  that  propagates  longitudinally  along  the  structure  at  about  35  m/s.  The 
hydrophones  nearest  to  the  module  ends  are  thus  less  than  one -half  wavelength 
away  from  the  ends  at  frequencies  less  than  40  Hz  while  at  the  same  frequency 
the  next  hydrophones  are  several  wavelengths  away.  This  structural  wave  is 
damped  by  the  hosewall  material  loss  tangent  and  is  significantly  damped 
within  a  couple  of  wavelengths.  Furthermore,  those  waves  are  efficiently  gen¬ 
erated  by  impedance  discontinuities  in  the  hosewall  such  as  occur  at  module 
end  terminations.  We  thus  surmise  that  such  structural  waves  are  generated  at 
the  acoustic  module  ends  and  are  observed  in  the  end  hydrophones  but  are  not 
seen  on  the  next  hydrophones  because  of  structural  damping.  This  cannot  be 
directly  measured  because  the  same  wave  must  be  measured  by  at  least  one 
hydrophone  pair.  There  were  other  times  when  all  of  the  hydrophones  in  the 
array  appeared  to  have  significant  low  frequency  signals  that  were  presumably 
vibration  induced.  Unfortunately,  when  this  condition  was  present  the  signals 
were  large  enough  that  an  electronic  overload  in  the  IPV  occurred  and  defini¬ 
tive  analysis  of  the  characteristic  of  the  vibration  could  not  be  made. 

Overload  protection  in  the  form  of  spectral  equalization  incorporated  in 
the  original  design  took  the  form  of  a  single-pole  high-pass  filter  section  in 
the  IPV  with  a  "corner"  fiequency  of  40  Hz.1  In  nearly  all  of  the  data 
obtained  from  this  deployment  it  was  seen  that  this  protection  was  inadequate. 
Figure  11  is  a  contour  plot  of  the  wavenumber- frequency  spectra  (equation  (2)) 
of  a  typical  experimental  trial  at  low  frequency.  The  wavenumber- frequency 
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Figure  11.  Aperture  7  Wavenumber- Frequency  Spectra  Illustrating 
Subsonic  Data  Contamination 


ACOUSTIC  LIMIT 


spectra  is  used  to  separate  the  sonic  and  subsonic  contributions  to  the 
hydrophone  outputs.  The  limits  of  the  "sonic  cone"  are  shown  as  lines  super¬ 
imposed  on  the  figure  using  the  velocity  at  the  face  of  the  array  from  figure 
8.  The  figure  shows  that  all  of  the  energy  above  25  Hz  is  acoustic.  Below  25 
Hz  significant  subsonic  energy  appears  although  the  sonic  levels  are  higher 
than  the  subsonic  levels.  The  character  of  the  subsonic  energy  is  further 
discussed  below. 

One  can  only  speculate  on  the  source  of  the  structural  vibration,  but  a 
likely  candidate  is  vortex/wake- induced  strumming  of  the  structure  caused  by 
ocean  currents.1  Using  the  ocean  current  values  inferred  from  figure  2,  one 
can  note  that  the  vibration  signals  were  infrequently  present  when  the  current 
was  less  than  1  to  3  knots.  A  particular  value  of  array  tilt  was  not  a  per¬ 
fect  predictor  of  the  existence  of  a  vibration- induced  signal  with  a  current 
less  than  1  to  2  knots.  However,  if  the  current  exceeded  1  to  2  knots,  then 
vibration  signals  were  always  present,  and  all  of  the  data  were  overloaded  if 
the  current  exceeded  2  to  3  knots. 

It  is  useful  in  the  design  of  a  vertical  array  to  know  whether  structural 
strumming  outside  the  actual  acoustic  aperture  is  contributing  nonacoustic 
signals.  In  this  discussion  "outside  the  acoustic  aperture"  refers  to  the 
structure  above  and  below  the  upper  and  lower  VIMs  shown  in  figure  1.  If  that 
portion  of  the  structure  is  generating  significant  vibration,  it  would  be  nec¬ 
essary  for  the  system  to  propagate  that  energy  as  a  longitudinal  extensional 
wave  because  transverse  vibrational  waves  in  water  are  rather  highly  damped  by 
the  fluid.  The  passage  of  longitudinal  waves  through  the  acoustic  structure 
would  be  observed  by  the  hydrophones  because  with  the  longitudinal  strength 
member  in  the  wall,  the  hosewall  would  couple  a  pressure  wave  to  the  hydroph¬ 
ones  by  the  hosewall  Poisson  effects.1  This  is  probably  the  source  of  the  570 
m/s  wave  (marked  on  figure  11)  observed  at  17  Hz,  which  is  the  frequency  for 
which  the  module  is  two  extensional  waves  long.  The  frequency  resonance  seen 
in  the  figure  is  thus  the  modal  response  of  the  structure  to  the  extensional 
wave.  The  existence  of  a  well-defined  wave  speed  was  not  typical,  but  the 
increase  in  subsonic  energy  below  20  to  25  Hz  was  always  the  case.  In  refer¬ 
ence  1  the  role  of  VIMs  as  a  low-pass  filter  to  longitudinal  vibration  was 
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discussed  and  the  "corner"  frequency  of  that  passband  was  calculated  to  be  23 
Hz.  A  better  design  would  have  had  the  extensional  two -wavelength  resonance 
above  the  VIM  cutoff.  The  conclusions  drawn  from  these  observations  are  that 
(1)  no  significant  nonacoustic  energy  is  present  in  the  measurement  band  of  40 
to  4000  Hz,  (2)  there  is  measurable  subsonic  energy  below  the  passband  of  the 
VIMs ,  i.e.,  25  Hz,  (3)  the  data  indicate  quite  significant  sonic  levels  at 
least  down  to  5  Hz,  and  (4)  the  single -pole  high-pass  filter  in  the  IPV  is  not 
adequate,  by  itself,  to  prevent  low  frequency  sonic  and  subsonic  energy  from 
occasionally  overloading  the  IPV  electronics.  Additional  equalization  is 
available  and  is  recommended. 

A  complete  discussion  of  the  preamplifier  noise  contribution  is  found  in 
reference  1.  It  was  anticipated  that  preamplifier  noise  could  be  a  limiting 
factor  for  wind  speeds  near  0  knots .  In  the  present  data  the  lowest  wind 
speed  case  was  4  knots,  which  produced  recording  levels  approximately  6  dB 
above  the  preamplifier  noise  contribution.  The  beamformer  output  data  at 
angles  below  the  horizontal  are  shown,  in  this  report,  to  have  relatively  low 
acoustic  sound  pressure  levels  as  a  consequence  of  the  source  being  near  the 
surface  which  required  a  bottom  reflection  to  arrive  at  angles  below  the  hori¬ 
zontal.  Figure  12  illustrates  the  beamformer  output  steered  "straight  down" 
for  the  lowest  wind  speed  data  recorded,  i.e.,  4  knots.  Superimposed  on  the 
measured  data  are  calculations  of  the  preamplifier  noise  contribution  to  the 
beamformer  output.  As  pointed  out  in  reference  1,  these  values  can  be  consid¬ 
ered  an  upper  bound  as  the  measured  points  typically  fell  about  3  dB  lower 
than  the  calculated  values.  A  measured  point  is  also  shown  in  the  figure. 
Comparison  of  these  two  curves  at  high  frequency  shows  that  preamplifier  noise 
is  contaminating  the  data. 
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Figure  12.  Lowes*-  Beamformer  Output  Levels  Measured  Compared 
With  Calculatec  Preamplifier  Noise  at  Beamformer  Output 
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4 . 0  EXPERIMENT  RESULTS 

4.1  TOTAL  SOUND  PRESSURE  LEVEL 

While  the  intent  of  this  experiment  was  to  measure  the  directional  spec¬ 
tra,  the  omnidirectional  sound  pressure  level  was  found  and  compared  with  his¬ 
torical  results.  The  hydrophones  used  had  no  significant  directional  pattern 
at  the  highest  frequency  of  the  measurements,  so  the  major  diagonal  terms  of 
the  csd  matrix  were  the  autospectral  density  of  the  total  field  sound  pressure 
level.  The  average  of  the  major  diagonal  terms  was  found  at  each  frequency 
of  interest.  The  results  for  a  selection  of  frequencies  are  shown  in  figure 
13,  displayed  as  functions  of  the  logarithm  of  wind  speed.  The  lowest 
frequency  shown  is  634  Hz  because  the  data  at  lower  frequencies  than  that  were 
frequently  contaminated  by  acoustic  sources,  to  be  discussed,  that  are  not 
related  to  wind  speed.  Superimposed  on  this  plot  are  data  points  taken  from 
the  Caron- Sheffield  memorandum18  which  is  an  extensive  data  base  of  total 
acoustic  sound  pressure  levels  for  the  area.  As  expected,  the  agreement  with 
the  Caron-Sheffield  data  is  good. 

In  figure  14  the  hydrophone  sound  pressure  level  data  are  replotted  as 
functions  of  nondimens ional  friction  velocity,  i.e.,  surface  friction  velocity 
nondimens ionalized  by  the  minimum  phase  speed  of  the  cap illary -gravity  wave. 
This  approach  is  identical  to  that  discussed  by  Kerman,  1984,  and  the  algo¬ 
rithms  used  here  are  the  same  as  found  in  reference  19.  The  traditional 
wisdom  is  that  a  value  of  1  for  this  nondimens ional  variable  represents  inci¬ 
pient  whitecapping.  Each  of  the  data  sets  is  partitioned  into  "whitecaps  pre¬ 
sent"  (wp)  and  "whitecaps  not  present"  (wnp)  according  to  whether  the  nondi- 
mensional  friction  velocity  is  greater  or  less  than  1,  respectively.  A 
regression  line  is  fit  to  each  partition.  The  data  base  is  quite  limited  for 
the  wp  case  and  no  power  ?aw  exponent  can  be  reliably  reported.  Qualitatively 
one  can  only  report  that  the  wp  case  appears  to  have  a  much  lower  dependence 
on  the  friction  velocity,  as  has  previously  been  reported.19  The  wnp  case  has 
a  larger  data  representation,  and  a  power  law  exponent  was  found  for  each 
frequency.  The  values  of  the  exponents  found  varied  randomly  from  1.4  to  2.2, 
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Figure  13.  Comparison  of  Autospectral  Density  With  Historical  Data 
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Figure  14.  Hydrophone  Autospectral  Density  Dependence  on 
Nondimens ional  Friction  Velocity 
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with  an  average  value  of  1.7.  This  is  lower  than  the  3.0  value  reported  by 
Kerman,19  for  reasons  to  be  discussed. 

4.2  DIRECTION - OF -ARRIVAL  SPECTRA  (DOAS) 

The  most  complete  view  of  the  results  of  a  particular  trial  is  the  eleva¬ 
tion  angle -frequency  DOAS.  The  frequency  range  of  40  to  4000  Hz  is  covered 
with  seven  nested  antennas.  Thus  the  multidimensional  spectra  consists  of  a 
composite  of  the  results  of  the  seven  apertures.  An  example  of  such  a  surface 
is  shown  in  figure  15.  Only  the  upper  octave  of  each  acoustic  aperture  is 
included  in  the  presentation.  The  figure  contains  both  a  surface  plot  of  the 
spectra  and  a  pseudo-colored  contour  plot  of  the  spectra.  Comprehending  the 
numerous  details  of  the  data  set  requires  comparing  both  views  of  the  three- 
dimensional  spectra  surface.  The  latter  is  frequency  "prewhitened"  to  reduce 
the  dynamic  range  of  the  plot.  The  "prewhitening"  is  accomplished  by  inte¬ 
grating  the  surface  with  respect  to  the  elevation  angle  at  each  frequency.  A 
third-order  regression  curve  is  fit  to  the  results  of  the  integral,  and  the 
resulting  polynomial  is  then  used  to  "whiten”  the  surface. 

A  spectra  similar  to  figure  15  is  available  for  each  of  the  25  data  sets. 
While  each  multidimensional  spectra  has  certain  unique  characteristics,  they 
all  fall  into  a  couple  of  categories,  examples  of  which  are  shown  in  figures 
15  through  18.  Figure  15  illustrates  three  characteristics  common  to  all 
data.  First,  the  six  peaks  occurring  in  the  downward- looking  direction,  i.e., 
180°,  are  the  aliasing  lotas  occurring  whenever  the  hydrophone  spacing  is 
exactly  one-half  wavelength.  Second,  the  "hump"  of  energy  centered  on  90° 
around  100  Hz  is  the  common-mode  noise  discussed  in  secc  m  3.  Third,  at  the 
frequencies  at  which  aliasing  first  occurs  there  is  a  discontinuity  in  the 
spectral  level  associated  with  the  change  in  aperture.  When  the  distribution 
of  energy  in  an  angular  direction  is  larger  than  the  spatial  resolution  of  the 
beamformer,  the  spectra  level  of  the  beamformer  will  increase  as  the  frequency 
decreases  because  the  beamwidth  of  the  aperture  increases.  At  the  demarcation 
between  apertures  the  beamwidth  changes  abruptly  from  a  four-wavelength 
antenna  to  a  two -wave length  antenna  with  a  subsequent  3  dB  drop  in  spectra 
level.  This  phenomena  is  better  observed  in  figure  19,  which  is  an  elevation 
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Figure  19.  DOAS  "Straight  Up"  Cut  Illustrating  a  Diffuse  Acoustic  Source 
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angle  "cut"  through  the  surface  at  0°,  i.e.,  "straight  up."  These  demarca¬ 
tions  are  not  present  if  the  spatial  extent  of  the  acoustic  source  is  less 
than  the  spatial  resolution  of  the  antenna.  This  is  illustrated  in  figure  20, 
which  is  an  elevation  angle  "cut"  at  82°  of  the  DOAS  shown  figure  20.  In 
figure  18  a  broadband  "target"  is  evident  at  that  elevation  angle.  Because 
the  "target"  is  spatially  compact,  the  aperture  demarcations  are  not  present 
in  figure  20.  These  three  characteristics  will  be  evident  to  some  degree  in 
all  of  the  figures  and  are  a  detraction  from  the  important  features  of  the 
figures. 

The  important  feature  of  figure  15,  which  is  characteristic  of  all  except 
the  very  low  wind  speed  data,  are  the  levels  above  the  horizontal,  i.e.,  90°, 
are  significantly  higher  than  below  it,  with  a  broad  peaking  of  the  spectra  at 
an  elevation  angle  of  about  40°.  This  is  characteristic  of  a  plane  or  a  thin 
layer  of  sources,  such  as  are  present  near  the  sea  surface.  The  beamformer 
coverage  of  this  acoustic  source  area  increases  as  the  elevation  angle  becomes 
larger,  i.e.,  further  from  the  zenith.  However,  the  received  source  level  is 
reduced  by  the  larger  source- to-receiver  distance  as  the  elevation  angle 
becomes  larger.  These  two  opposing  phenomena  cause  the  maximum  in  the  DOAS. 

In  calculating  the  directional  spectra,  the  received  acoustic  levels  are  nor¬ 
malized  by  the  solid- angle  coverage  of  the  beamformer,  which  removes  this 
peak.  Figure  16  is  quite  similar  to  figure  15  except  that  significantly 
larger  values  of  spatially  diffuse  energy  are  present  below  100  Hz.  This  pre¬ 
sumably  local  shipping-  and  industrial-source  acoustic  radiation  is  quite 
frequently  observed  contaminating  much  of  the  data  taken  below  100  Hz.  Also 
note  in  this  figure  a  narrowband  target  at  80°  elevation  angle  and  3500  Hz. 
Figure  17  is  an  example  of  a  low  wind  speed  data  set,  i.e.,  no  whitecaps  pre¬ 
sent.  The  general  pattern  of  the  DOAS  is  not  changed,  only  the  levels  are 
lower.  However,  note  that  at  high  frequency  and  at  an  elevation  angle  of  80° 
there  is  a  broadband  "target."  In  figure  18  this  "target"  is  better  defined. 
Figure  18  displays  the  lowest  level  data  obtained.  Here  the  aliasing  lobes 
have  completely  disappeared  because  preamplifier  noise  that  is  contaminating 
the  diffuse-source  noise  has  the  same  value  at  0*  and  180°.  A  broadband 
discrete  source  is  clearly  evident  at  82°.  This  shallow  angle,  8°  above  the 
horizontal,  is  most  likely  a  distant  source  arriving  along  a  bottom  bounce 
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Figure  20.  Low  Wind  Speed  DOAS  at  82°  Elevation  Angle  Illustrating 

a  Discrete  Acoustic  Source 
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path.  This  energy  is  nearly  always  present  when  the  wind  speed  is  low  enough 
to  observe  it  and  probably  should  be  considered  as  part  of  the  acoustic 
ambient  at  this  location.  There  are  two  possible  explanations  for  this 
energy.  This  elevation  angle  corresponds  almost  exactly  with  the  direction 
from  which  sounds  of  biological,  industrial,  or  shallow-water  wave  action  ori¬ 
gin  coming  from  the  top  of  the  escarpment  would  reach  the  measurement  site. 

An  alternative  source  is  acoustic  energy  in  the  TOTO  basin  that  is  efficiently 
coupling  to  the  measurement  antenna  by  the  waveguide  selection  process  of  the 
propagation.  The  TOTO  long-distance  propagation  is  characterized  by  a  bottom- 
surface  reflected  path  which  efficiently  propagates  at  only  a  narrow  selection 
of  angles.  The  permanent  refractive  character  of  the  media  limits  the  small 
grazing  angle  bottom  reflections,  and  the  bottom  reflection  loss,  which 
increases  with  increasing  grazing  angle,  limits  the  high-angle  arrivals  to 
relatively  short  range.  Thus  the  channel  selectively  limits  the  propagation 
mode,  i.e. ,  angles,  to  a  relatively  narrow  angular  window.  A  calculation  was 
made  of  the  received  vertical  angle  of  all  eigenrays  between  a  specific  source 
and  a  receiver  in  TOTO.  Figure  21  shows  the  received  angle  of  all  multipaths, 
i.e.,  eigenrays,  between  a  near-surface  source  and  a  213-meter-deep  receiver 
separated  by  38  km.  The  peak  of  the  structure  is  7°  from  the  horizontal, 
which  closely  matches  the  8°  shown  in  figures  17  and  18.  The  lack  of  symmetry 
about  the  horizontal  in  the  current  data  relative  to  figure  22  is  presumably  a 
result  of  a  multipath  being  "stripped"  by  the  sloping  approach  to  the  measure¬ 
ment  site.  Additional  support  for  this  argument  is  seen  in  figure  17,  which 
shows  an  easily  identifieo  narrowband  target  at  82°  and  3500  Hz.  This  is  a 
known  source  located  some  30  to  40  miles  from  the  measurement  location.  This 
identifies  the  preferred  propagation  path  as  being  identical  to  the  above 
argument.  The  frequency  spectra  of  the  beamformer  steered  to  the  energy  being 
discussed,  figure  23,  shows  a  frequency  spectral  shape  quite  similar  to  that 
of  „he  familiar  "snapping  shrimp"  described  by  Albers.20  Albers  documented 
sound  pressure  levels  in  the  vicinity  of  the  shrimp  "beds"  of  68  dB//ipa  in  a 
50  Hz  band,  although  there  is  some  ambiguity  in  the  description  with  regard  to 
the  measurement  bandwidth  This  means  a  spectral  level  in  excess  of  65 
dB//jiPa2/Hz  relative  to  our  measured  value  of  44  dB//pPa2/Hz  at  4000  Hz.  This 
level  difference  implies  that  the  source  is  local  because  the  propagation  loss 
in  a  typical  18  km  bottom,  surface/bottora  propagation  cycle  in  the  2  kHz  to 
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4  kHz  range  would  be  expected  to  be  larger  than  this  difference.  The  data  are 
thus  unable  to  choose  between  the  two  alternatives . 

For  either  explanation  we  conclude  that  the  acoustic  energy  arriving  via 
this  path  forms  a  nearly  stationary  contribution.  Figure  23  shows  the  DOAS  at 
82°  for  those  cases  in  which  the  presumed  "channel"  radiation  nearly  dominates 
the  frequency  range  shown.  For  the  times  in  which  measurement  was  possible, 
the  level  changed  by  less  than  4  dB.  The  fact  that  some  discontinuity  is 
observed  at  the  aperture  changes,  as  noted  in  the  figure,  indicates  that  there 
is  probably  a  spatially  diffuse  contribution,  indicating  that  the  range  of 
values  is  even  less  than  4  dB.  Figure  22  is  a  bearing  "cut"  of  a  selection  of 
measurements  chosen  to  illustrate  a  range  of  wind  speeds.  As  seen  in  the  fig¬ 
ure,  the  acoustic  data  will  be  affected  by  this  source  of  energy  below  a  4  m/s 
wind  speed.  This  character  of  the  data  will  be  discussed  again  when  we  pre¬ 
sent  the  directional  spectra  data.  The  local  maxima  in  the  data  at  180°  is  a 
result  of  the  aliasing  lobes. 

4.3  EIGENVECTOR  DECOMPOSITION 

As  discussed  in  section  3,  there  are  two  products  of  the  eigenstructure 
decomposition,  i.e.,  eigenvalues  and  eigenbeams.  When  the  eigenvalues,  which 
are  positive  scalars  for  a  Hermetian  matrix,  are  found  they  are  sorted  by  mag¬ 
nitude.  The  convention  we  will  use  is  that  the  largest  eigenvalue  is  always 
numbered  1.  Following  reference  6,  the  eigenvalues  are  scaled  so  that  their 
Pythagorean  sum  equals  the  average  of  the  major  diagonal  terms  in  the  csd  mat¬ 
rix,  i.e.,  the  total  acoustic  sound  pressure  level.  This  approach  follows 
from  the  orthogonality  of  each  of  the  terms  in  the  eigenstructure  decomposi¬ 
tion.  Four  examples  of  the  eigenstructure  decomposition  of  the  data  are  dis¬ 
played  in  figure  24.  Each  block  in  that  figure  shows  the  distribution  of  the 
acoustic  autospectral  density  by  eigenvector  for  seven  frequencies.  Pecall 
that  each  eigenvector  represents  the  spectral  density  contribution  of  a  par¬ 
ticular  statistically  independent  component  of  the  vertical  ambient  pressure 
field.  There  are  three  characteristics  of  these  plots  that  should  be  noted. 
First,  and  most  important,  there  is  an  underlying  regularity  to  the  distribu¬ 
tion.  Second,  the  regularity  is  broken  in  several  cases  by  the  appearance  of 
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Figure  24.  Examples  of  Eigenvalue  Decomposition 
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an  anomalously  large  value  in  the  first  or  second  eigenvalue.  When  the  eigen- 
beams  of  these  particular  eigenvalues  are  examined  we  find  that  they  are  asso¬ 
ciated  with  acoustic  plane  waves  or  contamination,  such  as  common-mode  noise, 
which  masks  as  a  plane  wave.  These  two  characteristics  illustrate  how  the 
eigenvector  decomposition  partitions  the  acoustic  field  into  discrete  sources 
(plane  waves)  having  a  discontinuous  eigenvalue  spectra  and  diffuse  sources  (a 
spatially  uniform  distribution  of  energy)  which  have  a  continuous  eigenvalue 
spectra.  The  third  characteristic  of  the  figures  that  needs  to  be  pointed  out 
is  that  the  79.81  and  3177.31  Hz  curves  "fall  off"  faster  with  increasing  eig¬ 
envalue  number  than  the  other  frequencies.  This  is  caused  by  the  reduction  in 
the  number  of  independent  eigenvectors  caused  by  a  reduced  nondimens ional 
acoustic  aperture.  The  number  of  independent  plane  waves  that  can  be  resolved 
by  an  antenna  is  twice  the  nondimens ional  aperture,  i.e.,  the  ratio  of  the 
physical  aperture  divided  by  the  wavelength.  At  80  Hz  the  aperture  was 
reduced,  relative  to  the  other  apertures,  because  of  the  loss  of  the  two  end 
hydrophones  caused  by  the  hydroacoustic- induced  voltage  overload  discussed  in 
section  3.  The  80  Hz  aperture  was  further  reduced  by  the  failure  of  hydro¬ 
phone  number  31,  so  this  low  frequency  aperture  only  had  five  hydrophones. 

Thus  there  are  only  five  independent  eigenvectors  present.  The  same  holds 
true  for  the  3200  Hz  frequency,  only  for  a  different  reason.  In  this  case  the 
aperture  reduction  is  because  the  frequency  is  further  below  the  design 
frequency  of  that  aperture  than  is  true  at  the  other  frequencies.  This  can  be 
seen  by  comparing  each  analysis  frequency  with  the  design  frequency  of  the 
antenna  used  (table  1).  A.t  3200  Hz  the  nondimens ional  aperture  is  2.3  com¬ 
pared  with  3.6  for  the  other  apertures,  except  for  80  Hz  which  is  2.5. 

In  figure  24(a)  an  obvious  discrete  source  is  present  at  six  of  the  seven 
frequencies.  This  corresponds  to  the  D0AS  pictured  in  figure  18  which  shows 
that  the  low  frequency  plane  wave  is  the  common-mode  noise  and  at  high 
frequency  the  discrete  source  appears  as  a  broad  bandwidth  "target"  at  an  80° 
elevation  angle.  In  figures  24(b)  and  24(c)  the  only  two  discrete  sources  are 
due  to  the  common-mode  noise.  In  figure  24(d)  the  common-mode  noise  appears 
only  at  159.24  Hz  because  an  energetic  low  frequency  diffuse  source  appears  to 
be  present  at  80  Hz. 
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Next  we  describe  the  results  of  the  eigenbeam  analysis,  where  the  eigen- 
beams  are  calculated  according  to  equation  (6)  of  section  2.  Eigenbeams 
resulting  from  a  purely  diffuse  acoustic  field  occurred  in  two  distinctly  dif¬ 
ferent  forms,  as  illustrated  in  figure  25,  which  shows  a  three-dimensional 
function  with  both  elevation  angle  and  eigenbeam  number  as  independent  vari¬ 
ables.  The  dependent  variable,  i.e.,  the  eigenbeam,  is  an  arbitrary  scale 
expressed  in  dB.  As  before,  the  eigenvectors  and  thus  the  eigenbeams  are 
sorted  according  to  the  magnitude  of  the  eigenvector,  so  the  values  of  the 
peaks  decrease  monotonically  from  back  to  front  of  the  plots.  Figure  25(a) 
shows  a  pair  of  diverging,  back- to -front  ridges.  Each  eigenbeam  appears  as  a 
pair  of  mainlobes  except  the  first  one  in  which  the  pair  has  merged  into  a 
single  lobe.  Because  each  eigenbeam  represents  a  component  which  is  statisti¬ 
cally  independent  of  all  other  eigenbeams,  the  pattern  is  consistent  with  a 
distributed  sea  surface  source  and  its  bottom  reflection.  As  the  sources 
become  more  distant  from  the  antenna,  both  the  direct  and  bottom  bounce  occur 
at  shallower  angles  until  the  two  merge  into  a  single  arrival,  slightly  shal¬ 
lower  than  90° .  This  pattern  is  typical  of  the  data  at  low  frequencies 
although  it  was  also  observed  at  a  high  frequency  and  low  wind  speed  condi¬ 
tion.  Of  particular  interest  in  this  pattern  is  that  the  eigenbeam  associated 
with  the  near-horizontal  elevation  angle  has  the  largest  amount  of  energy.  If 
the  sources  are  monopole  in  nature  then  the  energy  in  the  near -horizontal 
beams  would  be  larger  than  at  the  steeper  angles  because  of  the  larger  source 
area  contributing  to  those  solid  angles  as  long  as  there  is  sufficient  source 
energy  in  that  direction,  as  would  be  the  case  for  a  monopole  source.21  The 
second  pattern  is  seen  in  figure  25(b).  Here  the  acoustic  field  decomposes 
into  eight  statistically  independent  directions  with  the  levels  monotonically 
decreasing  from  "straight  up"  to  "straight  down,"  indicating  no  dependence 
between  the  direct  arrivals  from  the  surface  and  (presumably)  the  reflections 
off  the  bottom.  Note,  however,  that  now  the  maximum  eigenvalues  are  associ¬ 
ated  with  the  energy  coming  from  directly  overhead  and  diminish  monotonically 
as  the  elevation  angle  goes  to  90°,  as  would  be  the  case  if  the  diffuse  source 
were  dipoles.  The  continual  decrease  in  eigenvalue  level  as  the  angle  contin¬ 
ues  through  180°  is  presumably  caused  by  bottom  reflection  losses.  This 
second  pattern  is  typically  observed  above  a  couple  hundred  hertz,  but  it  was 
occasionally  seen  at  low  frequency  and  high  wind  speeds.  To  the  extent  that  a 
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general  statement  can  be  made,  the  low  frequency  eigenbeams  were  like  figure 
25(a)  and  a  conversion  to  the  figure  25(b)  pattern  occurred  above  500  Hz. 

As  we  saw  in  figure  24,  a  plane -wave -like  component  was  identified  in  the 
eigenvalues  as  an  anomalous  level  for  the  first  or  second  eigenvalue.  Figure 
26  illustrates  the  eigenbeams  associated  with  the  plane -wave -like  components. 
There  are  two  basic  forms  that  these  components  take  in  the  eigenbeam  plots. 
The  first  kind,  shown  in  figure  26(a),  is  when  the  plane  wave,  i.e.,  maximum 
eigenvalue,  occurs  in  the  same  direction  that  the  natural  ambient  maximum 
would  occur.  Thus  the  general  pattern  of  the  picture  is  an  eigenbeam  distri¬ 
bution  resulting  from  a  diffuse  monopole  source  with  an  anomalously  large 
value  at  90°,  which  in  this  case  is  the  common-mode  noise.  In  figure  26(b) 
the  plane-wave-like  eigenbeam  does  not  "fall  in  line."  The  "natural"  maximum, 
i.e.,  f^r^t,  eigenbeam  is  0°,  not  82°  as  seen  in  the  figure.  The  energy  of 
the  eigenbeam  at  82°  contains  two  independent  acoustic  sources:  (1)  sea 
surface  contribution  from  that  angle,  and  (2)  a  plane  wave  from  that  angle, 
which  in  this  case  is  the  TOTO  "channel"  radiation  discussed  previously.  It 
Is  the  contribution  from  the  second  source  that  caused  the  anomalously  large 
eigenvalue.  There  are  two  important  points  to  be  made  from  this  analysis. 
First,  in  the  absence  of  a  plane  wave  acoustic  source  the  magnitude -ranked 
eigenvalues  form  a  smooth  curve  as  noted  in  reference  6.  The  corollary  to 
this  is  that  discontinuities  in  the  magnitude -ranked  e^benvalue  distribution 
mark  the  presence  of  an  independent  plane  wave  in  the  measured  acoustic  field. 
Second,  the  eigenbeam  plot  identifies  the  eigenvector  that  has  been  perturbed 
by  the  plane  wave.  These  two  points  combine  to  form  a  simple  procedure  for 
removing  the  plane  wave  from  the  csd  matrix.  This  procedure  will  be  used  in 
the  next  section. 

4.4  VERTICAL  DIRECTIONAL  SPECTRA 


A  first-order  approximation  of  the  directional  spectra  was  calculated 
from  the  DOAS  as  discussed  in  section  2.  Figure  27  is  an  example  of  this  cal¬ 
culation  done  for  the  narrowband  composite  DOAS  shown  in  figure  15.  One  can 
see,  by  comparing  with  figure  15,  that  this  is  a  better  measure  of  the  acous¬ 
tic  field  for  two  reasons.  First,  because  the  calculation  accounts  for  the 
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Figure  27.  Composite  Narrowband  Elevation  Angle  -  Frequency  Directional 

Spectra,  11.83  m/s  Wind  Speed 


TR  8579 


change  in  solid-angle  coverage  of  the  beamformer,  the  discontinuities  at  the 
aperture  changes  are  no  longer  present,  although  the  aliasing  errors  remain. 
Second,  for  the  same  reason,  the  broad  maxima  near  40°  elevation  angle  is 
moved  to  0°,  which  represents  the  directional  spectra  for  the  anticipated 
dipole-like  sea  surface  sound  source  structure. 

As  previously  discussed,  the  frequency  variation  of  the  spectra  is  rela¬ 
tively  smooth  and  thus  there  are  advantages  to  synthesizing  proportional  band¬ 
width  spectral  estimates.  Note  that  the  dependent  variable  remains  a  spectral 
density  function,  not  band  levels,  but  the  spectral  estimates  are  now  deter¬ 
mined  by  all  of  the  energy  in  the  proportional  bandwidth,  which  in  the  present 
case  will  be  one- third-octave  bands.  This  process  significantly  reduces  the 
amount  of  data  present  in  the  data  base.  The  uniformity  of  the  data  within  a 
one- third- octave  bandwidth  with  frequency  was  quite  good  above  200  Hz.  It  was 
not  always  as  good  below  200  Hz  because  the  data  often  exhibited  "tonal"  com¬ 
ponents.  In  these  cases  the  basic  assumption  is  compromised.  We  will  con¬ 
tinue  to  draw  attention  to  this  fact  as  it  affects  our  interpretation  of  the 
data.  However,  the  rest  of  the  data  shown  in  this  section  employed  the  pro¬ 
portional  bandwidth  analysis. 

We  begin  our  examination  of  the  directional  spectra  by  viewing  its  wind 
speed  dependence.  This  is  done  in  the  same  manner  that  the  wind  speed  depen¬ 
dence  of  the  total  sound  pressure  level  was  examined  in  section  4.1,  i.e.,  by 
fitting  a  power  law  model  to  the  nondiraensional  friction  velocity  separately 
for  the  "with"  and  "without"  whitecaps  present  cases.  However,  in  this  case 
the  dependent  variable  is  the  directional  spectra,  not  the  total  sound  pres¬ 
sure  level.  This  process  was  repeated  at  eight  logarithmically  spaced 
frequencies  and  nine  elevation  angles.  Figure  28  illustrates  examples  of  this 
description.  Common  to  all  of  the  analysis  was  the  result  that  the  direc¬ 
tional  spectra  was  independent  of  the  nondimensional  friction  velocity  when 
whitecaps  were  present.  Because  of  the  meager  amount  of  data  for  this  case, 
the  conclusion  should  not  be  considered  quantitatively  correct.  However,  the 
result  is  qualitatively  consistent  with  the  results  found  in  reference  19  for 
the  total  sound  pressure  level  and  the  results  of  reference  22  for  the 
inferred  directional  spectra.  When  whitecaps  were  not  present,  the  power  law 
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Figure  28.  Examples  of  the  Wind  Speed  Dependence  of  the  Directional  Spectra 
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is  seen  to  vary  significantly  with  frequency  and  elevation  angle.  By  tabulat¬ 
ing  the  power  law  fit  to  the  data  as  a  function  of  both  elevation  angle  and 
frequency,  one  can  view  the  relative  wind  speed  dependence  of  the  directional 
spectra  as  a  function  of  frequency  and  elevation  angle  for  the  nonwhitecapping 
case.  This  construction  is  shown  in  figure  29.  The  system  common-mode  noise 
contamination  is  removed  by  substituting  for  the  measured  values  at  90°  eleva¬ 
tion  angle  in  the  frequency  range  of  50  to  600  Hz  with  values  consistent  with 
the  values  in  the  vicinity  of  the  contaminated  region.  There  are  two  obvious 
features  in  the  figure.  First,  the  data  above  and  below  300  Hz  show  signifi¬ 
cantly  different  wind- speed  dependence.  This  is  an  important  point  that  will 
be  returned  to.  In  the  next  paragraph  we  will  discuss  evidence  that  prior  to 
the  onset  of  whitecapping  there  appears  to  be  a  lack  of  wind- related  noise 
sources  below  300  Hz  relative  to  above  300  Hz.  The  complex  spatial  and 
frequency  content  of  the  acoustic  field  below  300  Hz  is  evident  in  figures  16 
through  18.  Little  more  can  be  said  than  that  the  acoustic  ambient  below  300 
Hz  is  importantly  influenced  by  local  human  activity.  The  second  dominant 
feature  of  figure  29  is  the  asymmetrical  "saddle"  shape  in  the  elevation  angle 
dependence  above  300  Hz.  The  high  side  of  the  "saddle"  is  the  zenith  of  the 
array.  The  power  law  exponent  of  2.8  to  3  is  consistent  with  past  experi¬ 
ence.19,22  The  other  side  of  the  "saddle,"  i.e.,  looking  down  at  bottom 
reflections,  has  consistently  smaller  exponents,  i.e.,  1.5  to  2.0.  This  indi¬ 
cates  that  the  pressure  field  coming  from  that  direction  is  not  simply  the 
reflection  of  the  sea  surface  sound  but  includes  a  non- sea  surface  generated 
component  such  as  biological  activity.  The  lowest  wind  speed  dependence  in 
the  directional  spectra  above  300  Hz  is  seen  to  exist  from  approximately  70® 
to  100® .  The  likely  explanation  for  this  is  the  combination  of  escarpment 
radiation  and  "channel"  radiation  discussed  in  section  4.2.  In  summary,  the 
acoustic  ambient  at  all  frequencies  below  approximately  300  Hz  is  importantly 
influenced  by  local  human  activity  and  low  wind  speed  dependence  without 
whitecaps  present,  and  above  300  Hz  the  energy  coming  from  the  sea  surface  is 
dominated  by  a  cubic  wind  speed-dependent  mechanism.  Near  the  horiz  atal 
direction  and  in  downward  looking  directions  the  wind  speed  dependence  is 
reduced  by  the  bottom  reflection  and  there  is  a  clear  non-wind  speed  component 
present . 
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Note  that  the  calculation  of  the  directional  spectra  performed  in  this 
study  assumes  that  the  directional  spectra  varies  slowly  relative  to  the  beam-  1 

former  pattern  function  mainlobe  width.  This  would  not  be  the  case  when  a 
strong  plane  wave  is  present  in  the  measured  pressure  field.  Thus  all  of  the 
directional  spectra  presented  in  this  section  have  been  processed,  as  dis¬ 
cussed  in  section  4.3,  to  remove  plane  wave  "contamination."  This  includes  1 

the  "common-mode  noise,"  discussed  in  section  3,  which  appears  in  the  data  as 
a  plane  wave  at  90°.  Directional  spectra  with  all  plane  waves  removed  will  be 
labeled  as  "corrected  directional  spectra."  Displaying  the  data  is  difficult 
because  the  space  of  interest  has  four  dimensions,  i.e.,  spectral  level,  < 

frequency,  elevation  angle,  and  wind  speed.  Our  approach  will  be  to  show  the 
three-dimensional  frequency- elevation  angle  directional  spectra  with  wind 
speed  as  a  parameter.  Five  examples  of  the  data  were  chosen  to  represent  the 
range  of  wind  speed  data  observed.  Topographical  and  surface  plots  of  the  1 

figures  are  shown  in  figures  30  through  34.  Note  that  the  topographical  plots 
are  "whitened"  in  the  frequency  dimension  to  reduce  the  dynamic  range  but  the 
surface  plots  are  not  prewhitened.  There  are  several  features  of  this 

sequence  of  figures  worth  discussing.  As  the  wind  speed  increases,  the  energy  i 

at  frequencies  above  300  Hz  increases  more  rapidly  than  the  energy  below  300 
Hz.  Prior  to  incipient  whitecapping,  i.e.,  about  6.5  m/s,  the  frequency- 
elevation  angle  space  contains  two  types  of  regions.  One  region  has  the  maxi¬ 
mum  value  indicated  as  a  ridge  that  peaks  at  angles  nearly  the  same  as  in  the  { 

case  prior  to  the  removal  of  the  plane  wave.  The  result  acts  as  though  the 
removal  was  incomplete.  The  reason  for  this  is  that  the  component  still  pre¬ 
sent  after  the  plane  wave  removal  appears  from  the  beam  response  and  the  eig¬ 
envector  decomposition  to  be  spatially  diffuse.  This  is  not  a  surprising  con-  I 

elusion  for  either  the  "escarpment"  radiation  or  the  "channel"  radiation 
sources.  However,  it  should  not  be  the  case  for  a  "common-mode"  noise.  Thus 
we  conclude  that  a  spatially  diffuse  acoustic  source  is  present  which  at  times 
has  frequency  high  spots  that  are  harmonically  related.  It  was  coincidental  I 

that  the  "common  mode"  noise  peaked  on  this  low  frequency  ridge.  As  the  wind 
speed  increases,  the  nearly  symmetrical  pattern  previously  described  is 
replaced  with  a  nonsymmetL ical  pattern  showing  the  dominance  of  the  surface - 
generated  acoustic  radiation.  A  striking  feature  of  figures  31  and  32  is  the  i 

movement  of  one  pattern  to  the  other  near  300  Hz.  Finally,  we  note  that  when 
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whitecapping  is  present  the  asymmetrical  pattern  dominates  nearly  the  entire 
two  decades  of  frequency,  indicating  the  efficient  sound  sources  present  in 
the  whitecapping  process.  Thus  below  300  Hz  the  sea  surface  radiation  process 
is  quite  different  than  it  is  above  300  Hz  when  whitecaps  are  not  present. 
However,  with  whitecapping  the  mechanisms  below  and  above  300  Hz  are  very  sim¬ 
ilar  . 


Figures  35  through  38  illustrate  elevation  angle  "cuts"  of  the  previous 
three-dimensional  surfaces  of  the  corrected  directional  spectra  for  selected 
frequencies.  Because  the  frequencies  chosen  were  always  near  the  frequency 
for  which  the  hydrophone  spacing  was  one -half  wavelength,  the  resulting  spec¬ 
tra  is  biased  by  the  aliasing  error  of  the  array.  For  this  reason  the  biased 
values,  which  are  in  the  vicinity  of  180°,  have  been  removed.  Data  represent¬ 
ing  different  wind  speeds  are  overplotted  in  each  graph.  These  graphs  illumi¬ 
nate  the  discussions  of  the  previous  paragraph.  Note  in  particular  the  sig¬ 
nificant  change  in  the  downward-propagating  energy  when  whitecaps  were  pre¬ 
sent,  i,e.,  wind  speed  greater  than  6.5  m/s,  at  frequencies  less  than  300  Hz 
relative  to  greater  than  300  Hz. 

4.5  ANISOTROPIC  GAIN 

A  primary  function  of  hydrophone  arrays  is  the  detection  and  measurement 
of  underwater  acoustic  radiation.  A  measure  of  the  spatial  gain  that  one 
achieves  with  this  structure  is  found  by  taking  the  ratio  of  the  total  acous¬ 
tic  sound  pressure  level  to  the  beamformer  output  calibrated  to  properly 
measure  the  sound  pressure  level  of  a  plane  wave.  This  ratio,  for  a  properly 
designed  structure,  is  of  the  order  of  the  number  of  hydrophores  at  the 
frequency  for  which  the  hydrophones  are  spaced  an  acoustic  one -half  wavelength 
apart.  In  the  present  case  the  gain  is  between  8  to  9  dB,  depending  on  the 
spatial  weighting  function  used.  This  is  the  spatial  gain  of  the  system  in  an 
isotropic  ambient  field  and  is  referred  to  as  the  directivity  index.  The 
directivity  index  is  a  function  of  both  frequency  and  beamformer  steering 
angle;  however,  it  is  only  mildly  dependent  on  these  variables,  changing  by 
only  3  dB  over  the  180°  steering  range  and  an  octave  of  frequency  range. 
Clearly,  the  actual  spatial  gain  will  depend  sensitively  on  the  steered  direc- 
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tion  of  the  beamformer  in  an  anisotropic  acoustic  ambient.  In  this  section  of 
the  report  we  separate  the  spatial  gain  associated  with  the  anisotropic  pres¬ 
sure  field  by  normalizing  the  measured  spatial  gain  by  the  calculated  direc¬ 
tivity  index  at  each  frequency.  The  result  of  this  calculation  we  term  aniso¬ 
tropic  gain.  Such  a  presentation  amounts  to  a  transform  of  the  data  given  to 
another  observation  space,  i.e.,  directional  spectra  to  anisotropic  gain.  The 
independent  variables  remain  the  same.  Thus  the  data  shown  in  figures  30 
through  38  are  shown  as  anisotropic  gains  in  figures  39  through  47.  The  fig¬ 
ures  quantify  the  obvious;  i.e.,  high  spatial  gains  are  found  by  steering  to 
regions  of  low  acoustic  ambient  radiation.  At  low  wind  speeds  this  is 
achieved  by  steering  either  above  or  below  the  horizontal.  As  the  wind  speed 
increases  then,  depending  on  the  frequency,  gains  greater  tnan  the  directivity 
index  are  observed  from  the  horizontal  direction  on  down  to  180° .  Once  the 
wind  speed  becomes  high  enough  to  establish  this  asymmetrical  pattern,  then 
the  anisotropic  gain  becomes  relatively  insensitive  to  further  wind  speed 
increases.  The  increased  gain  obtained  by  steering  downward  is  quite  signifi¬ 
cant,  i.e.,  5  to  6  dB  over  the  9  dB  directivity  index.  Thus  the  spatial  gain 
of  an  8 -hydrophone  array  exhibits  the  gain  of  a  25 -hydrophone  array. 
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Figure  39  Corrected  Anisotropic  Gain  Versus  Frequency  and  Elevation 

Angle  for  2.06  m/s  Wind  Speed 
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Figure  40.  Corrected  Anisotropic  Gain  Versus  Frequency  and  Elevation 

Angle  for  2.83  m/s  Wind  Speed 
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Figure  42.  Corrected  Anisotropic  Gain  Versus  Frequency  and  Elevation 

Angle  for  7.7  m/s  Wind  Speed 
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Figure  43.  Corrected  Anisotropic  Gain  Versus  Frequency  and  Elevation 

Angle  for  11.83  m/s  Wind  Speed 
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SUMMARY 


This  report  documents  the  data  processing  techniques  and  the  results  of 
the  initial  deployment  of  a  wideband  (40  to  4000  Hz)  vertical  acoustic  antenna 
system  in  The  Tongue  of  the  Ocean,  The  Bahamas,  at  a  site  adjacent  to  the 
Andros  Island  eastern  escarpment.  The  instrumentation  system  consisted  of 
seven  octavely  nested  four-wavelength  acoustic  antennas  covering  the  above 
frequency  range  and  a  VCR-based  subsurface  recording  system  that  carried  out  a 
programmed  sequence  of  recordings.  The  purpose  of  the  experiment  was  the  ini¬ 
tial  field  test  of  a  deployable  acoustic  monitoring  system  (DAMS) .  The  suc¬ 
cess  of  the  experiment  made  available  a  data  set  that  allowed  the  examination 
of  the  vertical  spatial  structure  of  the  acoustic  ambient  of  the  location. 
While  there  were  aspects  of  the  data  unique  to  the  location,  the  bulk  of  the 
information  appears  to  apply  to  the  more -generic  problem  of  characterizing  the 
directional  spectra  resulting  from  a  sea  surface  sound  source  structure. 


The  two-month  experiment  period  of  September  and  October  1988  made  avail¬ 
able  25  data  sets  that  adequately  covered  the  "no  whitecaps  present"  case  but 
only  sparsely  covered  the  "whitecaps  present"  case.  The  measured  total  sound 
pressure  levels  compared  favorably  with  previous  experience,  indicating  that  a 
unity  value  of  the  surface  friction  velocity  nondimensionalized  by  the  minimum 
phase  velocity  of  the  surface  capillary-gravity  wave  field  marked  a  change  in 
the  acoustic  ambient  pressure  field  caused  by  incipient  whitecapping.  The 
power  law  fit  to  the  "no  whitecaps  present"  case  was  on  the  average  only  1.7 
relative  to  the  3.0  reported  in  reference  19.  However,  the  power  law  exponent 
fit  to  the  directional  spectra,  in  the  direction  of  the  surface,  above  300  Hz 
varied  from  2.8  to  3.1.  Upward- propagating  energy  had  exponents  that  varied 
from  1  to  2 ,  depending  on  the  elevation  angle  and  frequency.  Thus  even  at 
frequencies  above  300  Hz  the  sea  surface  sound  was  not  the  only  source  of 
acoustic  energy  present.  Below  300  Hz  the  power  law  exponent  was  always  less 
than  1.  There  appeared  to  be  two  reasons  for  this  lack  of  wind  speed  depen¬ 
dence.  First,  this  frequency  range  regularly  contained  narrowband  acoustic 
energy,  presumably  the  result  of  industrial  and  boating  activity  in  the  area. 
Underlying  these  randomly  occurring  events  was  a  spatially  diffuse  acoustic 
background  with  a  broad  maxima  in  the  horizontal  direction  that  exhibited 
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little  wind  speed  dependence.  Second,  prior  to  the  onset  of  whitecapping 
there  was  little  to  no  wind  speed  dependence  in  the  data.  Once  whitecapping  I 

occurred,  the  levels  jumped  noticeably.  This  appears  to  indicate  a  difference 
in  the  acoustic  sea  surface  source  mechanism  above  and  below  300  Hz  with  no 
whitecaps  present.  With  whitecaps  present,  the  acoustic  source  mechanism 
appeared  uniform  over  the  entire  two  frequency  decades  of  the  measurement,  I 

exhibiting  the  usual  broad  maxima  near  500  Hz,  which  has  been  attributed  to 
spray  impact  on  the  surface.23 

The  measured  vertical  directional  spectra  display  two  patterns  in  an  ele¬ 
vation  angle -frequency  space.  The  first  pattern  consists  of  a  ridge  of  energy 
at  or  near  the  horizontal  direction  which  drops  off,  i.e.,  less  energy,  mono- 
tonically  and  nearly  symmetrically  in  either  the  up  or  down  direction.  The 
entire  pattern  is  inversely  proportional  to  frequency.  The  second  pattern  is 
quite  nonsymmetrical  about  the  horizontal  with  a  broad  maxima  in  the  direction 
of  downward-propagating  energy  and  a  minima  in  the  direction  of  upward- 
propagating  energy.  The  latter  is  the  familiar  pattern  associated  with  a  sur¬ 
face  layer  of  independently  radiating  dipoles.  At  the  lowest  wind  speeds  of 
about  2  m/s  the  first  pattern  extends  over  the  entire  frequency  range.  As  the 
wind  speed  increases,  the  second  pattern  dominates  above  300  Hz  and  the  first 
pattern  dominates  below  300  Hz.  This  dual  pattern  continues  until  whitecap¬ 
ping  commences,  at  which  time  the  second  pattern  dominates  the  entire  two 
frequency  decades  of  the  measurement.  At  low  wind  speeds,  when  the  first  pat¬ 
tern  dominates  the  entire  frequency  range,  the  ridge  of  energy  above  500  Hz  is 
centered  at  8°  above  the  horizontal.  This  ridge  extends  to  4000  Hz.  There 
are  two  explanations  for  this  energy.  First,  this  is  the  direction  associated 
with  radiation  from  the  top  of  the  escarpment.  The  frequency  spectral  shape 
shows  increasing  levels  from  1500  to  4000  Hz,  which  is  consistent  with  "snap¬ 
ping  shrimp"  radiation.  This  explanation  indicates  a  local  sound  source. 

However,  it  is  known  that  the  Tongue  of  the  Ocean  acts  as  waveguide  that  effi¬ 
ciently  propagates  energy  via  a  bottom/surface  interacting  path  over  a  very 
narrow  angular  window.  Predictions  of  the  center  of  this  elevation  angle 
indicate  a  narrow  peak  at  about  7°  from  the  horizontal.  This  has  been  veri¬ 
fied  in  the  data  repeatedly  by  the  reception  of  sonar  pulses  propagated  from 
more  than  40  km  away  and  arriving  in  the  center  of  the  ridge  being  discussed. 
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This  explanation  indicates  that  long-distance  sound  sources  throughout  the 
#  "tongue"  would  communicate  with  the  measurement  on  this  ridge.  The  actual 

explanation  for  this  ridge  of  acoustic  energy  could  be  either  or  both  explana¬ 
tions  . 


The  basic  measurement  of  the  experiment  was  the  cross -spectral  density 
matrix  for  each  of  the  seven  vertical  hydrophone  arrays.  An  eigenvector 
decomposition  of  these  matrices  was  performed.  There  are  two  basic  results  of 
the  decomposition.  First,  the  eigenvalues,  sorted  by  magnitude,  indicate  a 
uniform  distribution  of  magnitudes  except  when  an  independently  occurring 
plane  wave  is  present  in  the  measurement.  Second,  the  eigenbeams,  when 
plotted  as  a  three-dimensional  surface  versus  elevation  angle  and  eigenvalue 
number,  form  two  patterns  consistent  with  the  two  directional  spectra  patterns 
discussed  above.  First  is  a  pattern  showing  a  diffuse  source  peaking  in  or 
near  the  horizontal  direction  and  uniformly  decreasing  in  both  the  up  and  down 
directions.  Each  eigenbeam  contains  two  peaks,  one  up  and  one  down.  A  pos¬ 
sible  interpretation  of  this  pattern  is  that  the  upward-  and  downward-propa¬ 
gating  energy  occurs  in  spatially  coherent  pairs.  The  second  pattern  indi¬ 
cates  a  spatially  diffuse  source  structure  with  a  maxima  for  downward-propa¬ 
gating  energy  and  a  minima  with  upward-propagating  energy.  The  uniform  dis¬ 
tribution  of  eigenvalues  was  used  to  remove  contaminating  plane  waves  by 
adjusting  the  level  of  the  eigenvalue  containing  the  plane  wave  to  be  consis¬ 
tent  with  the  uniform  distribution  of  eigenvalues  that  would  exist  in  the 
absence  of  the  plane  wave. 


Three-dimensional  surfaces  of  the  anisotropic  gain  were  plotted  as  a 
function  of  frequency  and  elevation  angle  with  wind  speed  as  a  parameter. 
Anisotropic  gain  is  defined  here  as  the  spatial  gain  of  the  vertical  hydro¬ 
phone  array  relative  to  the  gain  of  the  same  array  in  an  isotropic  noise 
field.  The  results  quantified  the  obvious;  i.e.,  beamformer  gains  steered  in 
the  downward  direction  can  realize  significant  improvements  over  the  directi¬ 
vity  index  of  that  array.  In  the  present  case  14  to  15  dB  of  spatial  gain  was 
measured  for  an  array  having  a  9  dB  directivity  index.  This  result  in  the 
downward  direction  was  relatively  independent  of  the  wind  speed. 
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